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O Abstract 

The Roth's two-pole approximation has been used by the present au- 
thors to study the effects of the hybridization in the superconducting 
^ ' properties of a strongly correlated electron system. The model used is 

I the extended Hubbard model which includes the d — p hybridization, the 

p-band and a narrow d-band. The present work is an extension of the 
previous Ref. [3]. Nevertheless, some important correlation functions nec- 
essary to estimate the Roth's band shift, are included together with the 
temperature T and the Coulomb interaction U to describe the supercon- 
' ductivity. The superconducting order parameter of a cuprate system, is 

obtained following Beenen and Edwards formalism. Here, we investigate 
' in detail the change of the order parameter associated to temperature, 

Coulomb interaction and Roth's band shift effects on superconductivity. 
The phase diagram with Tc versus the total occupation numbers nr, shows 
' the difference respect to the previous work. 
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' The high temperature superconductors discovered by Bednorz and MuUer 

^ I [1] are beUeved to be explained in the framework of the Hubbard model [2] since 

electron correlations are strong. 

In this work we have used the extended Hubbard model [3] with the Roth's 
method [4]. In the previous work [3], we improved Ref. [5] to also include su- 
perconducting properties following closely the approach introduced by Beenen 
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and Edwards [6]. Nevertheless, in Ref. [3], the Roth's band shift has been es- 
timated disregarding the temperature and the Coulomb interaction effects. In 
the present work, wc have included these effects in the Roth's band shift. Con- 
sequently, due to the U dependence, it is necessary to calculate more correlation 
functions than in the reference [3] . Therefore, we hope by using this procedure 
to get a more correct behavior of the superconductor order parameter (as a 
function of T) and the critical temperature T^. The Hamiltonian used is: 

i,a i 

i,o- i^J^o" 

+ Yt^t (dlpj^+pldj^) (1) 

where /z is the chemical potential. 

In order to study superconductivity by Roth's method, it is necessary to 

include "hole" operators in the set of "electron" operators (which describes the 
normal state) and evaluate anomalous correlation functions. As discussed in 
Ref. [3], considering the particular case with singlet pairing and the rf-wave 
symmetry, we get (di-a-di^) = and '^^{di^a-dia-) = 0, where the sum is over 
sites I which are nearest neighbors of i. Therefore, in the d-wave case [6], the 
superconducting gap is determined by the gap function. 



Ik = -Ik f 

^ k 



cos (/ca^a) - cos {kya)f -^j) f{Lo)Gll{ui)dui (2) 



where ,f{uj) is the Fermi function. The propagators G^^, with s = 1,2 or 3 
can be obtained as in Ref. [3]. Wc can readily obtain the correlation hmction 
nf„ = {d\^dicr) considering a nearest-neighbor model with t^j = for the z 
neighbors. In this way it is necessary only one value of n^^, which is nf^. The 
gap function given by equation (2), defined in Ref. [6], is solved sclf-consistently. 
The band shift Wka = + Wf, shifts the poles of the propagators G^^;^. In 
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this work we considered a k independent hybridization {V^ Y = {Vj^^Vj^' ) 
discussed in Ref. [5], where (...) is the average over the Brillouin zone and 
V^^{V^'^) are the Fourier transform of t^j{t^j). Consequently, the W^"^ term is 
k independent, and also, its temperature and Coulomb interaction dependence 
is quite small. Thus, our focus here is W^^ which is given by: 



1 + I- Ha 



1 - HaH-a 

(3) 
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Figure 1: as a function of the total occupation number nx (where nx = 
+ nl^). In (a) and (b), the dotted Hnes show the previous results from 
Rcf. [3] for U = 12|t''|. The solid lines show the behavior of Tc in the present 
approach. The figures (c) and (d) show the present results for U = 8|<''|. The 
units of the hybridization Vq"^ are electron-volts (eV). {f^ = —0.5eV.) 



where a and /3 are defined as: 

"- = T^^ and /3.>= ^.^(,_^.^) , (4) 

with nl^ = ng-o-- The correlation function rrija = {dl^n'j_„dja) is obtained 
from the propagator G^^ (see Ref. [3]) as: 



(5) 



The band shift given by equation (3) can be split into a k-dependent and a 
k-independent term. While the role of the k-independent term is only to shift 
the poles of the Green's function, the k-dependent part causes a flattening at 
the top of the lower band from around k = (tt, tt) until k = (tt, 0) as can be 
verified in Refs. [6, 8]. This flattening occurs because Wj^^ decreases when 
ef increases, as discussed in Ref. [6]. This numerical calculations agree with 
the Monte Carlo results obtained by Bulut et al [6, 8]. The flattening in the 
quasi-particle sub-bands (called Hubbard bands) leads to a band narrowing of 
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these sub-bands, increasing the gap originated by the Coulomb interaction U. 
Our numerical results show that the narrowing of the sub-bands decreases with 
increasing U. This fact is very important, because this leads to a change at the 
position of the chemical potential jj, which is relevant to obtain the behavior of 
the high-temperature superconductivity [8]. In the figures l(a)-l(b) are shown 
phase diagrams with the dotted lines corresponding to the results obtained in 
Ref. [3], where the effects of the temperature T and the Coulomb interaction 
are not included in W^^. In the figures l(c)-l(d) the value of the U is decreased. 
The solid lines show the results for Tc, with the effects of T and U included in 
the calculation of the IV^^ band shift. The results show a small increase of Tc 
when we consider in W^^ the effects of T and U. We believe that the increasing 
of Tc occurs as result of the change in the value of the chemical potential, 
as discussed above. The numerical results show also that in the interval of 
temperatures showed in figure 1 the effect of the temperature in W^^ is quite 
small. 
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